function [] = displayEstimates_acctCreation(runIds)

%%% Add paths
addpath('../general_funcs');

LL = 0;
NumParams = 0;

NumRuns = length(runIds);
for runidx = 1:NumRuns
	runId = runIds(runidx);
	
	%%%% Get input / output paths
	config = getConfig_acctCreation(runId);
	runFolder  = config.runFolder;

	% Load estimation results
	load(sprintf('%s/results.mat', runFolder), 'params_hat', 'params_ses', 'modelFits');
	
	%%% Load data to get paramNames
	load(sprintf('%s/runData.mat', runFolder), 'data');
	paramNames = get_param_names_acctCreation(data);
	NumObs = size(data.Y(:),1);
	clear data;
	
	%% Increment LL and NumParams
	LL = LL + modelFits.LL;
	NumParams = NumParams + length(paramNames);
	if runidx == 1;
		NumObs_all = NumObs;
	else
		if NumObs ~= NumObs_all
			error('Inconsistency in the number of observations across runs.')
		end
	end
	
	% Make table of estimates and display it
	M1 = reportSubmodelResults(paramNames, params_hat, params_ses, modelFits, {}, true);
	clear params_hat params_ses modelFits;
	disp(M1);
	
	% Write table to output file
	writetable(M1, sprintf('%s/results2.csv', runFolder), 'WriteRowNames', true);
end

AIC = 2*NumParams - 2*LL;
BIC = NumParams*log(NumObs) - 2*LL;

LL
AIC
BIC
NumParams
NumObs

M2 = table(LL, AIC, BIC, NumParams, NumObs);
writetable(M2, sprintf('%s/../modelFits.csv', runFolder), 'WriteRowNames', true);

end
